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The critical behavior of quarks undergoing phase transition to hadrons is consid- 
ered in the framework of the Ising model. It is found that spatial fluctuations 
do not alter the _F-scaling result obtained earlier in the Ginzberg-Landau formal- 
ism. For the study of chaos a new measure must be found for the quark-gluon 
system, for which the usual description appropriate for classical nonlinear systems 
is inapplicable. It is shown that the entropy indices fi q which characterize the 
fluctuations of spatial patterns, are as effective as the Lyapunov exponent in mea- 
suring chaos. When applied to the QCD parton showers, it is found that the quark 
jets are more chaotic than gluon jets. The analysis is highly appropriate for a 
quantitative treatment of the erratic fluctuations of multiparticle production in 
leptonic, hadronic and nuclear collision processes. It represents a step beyond the 
traditional intermittency analysis. 



1 Introduction 

The topics to be discussed here concern the nature of fluctuations of hadronic 
observables in high-energy collisions. The first is on quark-hadron phase phase 
transition, a subject that has been investigated from many approaches. Our 
emphasis is on the scaling behavior of what can be observed experimentally. 
Our earlier study of the subject has been on the basis of a mean-field theory. 
We now improve upon that by incorporating spatial fluctuations. 

The second topic is on the possibility of chaotic behaviors of a quark- 
gluon system. Since the number of degrees of freedom of such a system is 
not conserved, the traditional method of following the distance between two 
nearby trajectories of a classical system cannot be applied. Thus a major part 
of our effort in that problem is to identify an appropriate measure of chaos 
and to show that it is as good as the usual one when applied to the classical 
problems. With the new measure it is then possible to investigate the question 
of chaos for QCD parton showers. The significance of this line of investigation 
is in uncovering the fluctuating properties of multiparticle production that 
have hitherto been averaged over. The experimental determination of the 
key measures that quantify those properties will then present a challenge to 
theoretical models, which may not have been sufficiently accurate (at least in 
soft interaction) to reproduce them. 

The two topics will be discussed separately below, with more space given 
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to the latter, since it involves newer and more unfamiliar concepts. 
2 Critical Behavior 

In a series of papers 0i the consequences of quark-hadron phase transition 
in terms of hadronic observables have been studied in the framework of the 
Ginzburg-Landau formalism. It is found that the normalized factorial moments 
F q do not exhibit intermittency (i.e., the power-law behavior M Vq ) : but they 
do satisfy _F-scaling 

F q <xF$* , (1) 

where 

p q = (q-l) v , 1/ = 1.304 . (2) 

The value of the scaling exponent v, which bears no relationship to any of 
the usual critical exponents, is a definitive characteristic of second-order phase 
transition, in which the temperature is not a measurable quantity, and the 
observables are multiplicities of the quanta produced. Heavy-ion collision ex- 
periments have shown the validity of ([!]), but currently the value of v is not 
1.304, signifying the absence of quark-hadron phase transition. But in a laser 
experiment at the threshold of lasing, where the physics is known to be that of 
a second-order phase transition, it has been shown that (0) is valid to a high 
degree of accuracycl 

The Ginzburg-Landau theory is a mean-field theory, so the spatial proper- 
ties of the system are assumed to be smooth. To improve on that description, 
it is necessary to take into account the fluctuating character of hadronization. 
That has been done recently in two different directions. OneJs done by simula- 
tion on a 2D lattice using the Ising model for hadronizationa The other is also 
done on the lattice, but surface fluctuations (perpendicular to the lattice) are 
introduced with the thermal effects being constrained by the color confinement 
potential^ We discuss them in turn. 

2. 1 Spatial fluctuations simulated on the Ising lattice 

In the Ginzburg-Landau (GL) approach to the problem of phase transition we 
relate the order parameter <p(z) to the hadron density by p(z) — \<fr(z)\ . For 
a bin of area 5 2 we assume that p{z) is constant in the bin so the GL free 
energy is simplified to the form F[<f>] = 5 2 ia\4>\ 2 + b\(j)\ A \ . The multiplicity 
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distribution P n of hadrons for T < T c is then given 





4>{z)=A-^Y. s i > W 



where is the spin ±1 at site j, and z is now the location of the center of a 
cell of size A £ , which is < S 2 , but large enough to contain several sites. The 
hadron multiplicity at the ith cell is then 



where A is a scale factor that relates the quark density of the plasma at T c to 
the lattice site density in the Ising model. The bin multiplicity is n$ = n i; 
where N$ = (S/e) 2 . n 

By using the Wolff algorithm □ to simulate the bin configurations in the 
Ising model, we can calculate P n (S) and therefore F q (M), where M = (L/S) 2 , 
L being the lattice size. We have found that the strict scaling behavior, F q oc 
M Vq , is valid only at T = T c . In fact, T c can be determined by varying T until 
that behavior is manifested. However, for T <T C the more general ^-scaling 
behavior (1) is valid for a range of T. Furthermore, (2) is also valid, but the 
value of v depends on T. 

Our result is that v ~ 1.0 at T c but v becomes bigger at T < T c . For a 
range of T < T c , the values of v are between 1.0 and 1.6 so that the average 
value of v is about 1.3. This is a very satisfying solution to a dilemma posed 
by an earlier analytical result, where SatzQ found <p q = (q — l)/8 at T = T c in 
the 2D Ising model. It means that j3 q = f q /f2 = ? — 1, so that v — 1.0. Now, 
we see that the GL mean-field result of v = 1.304 is an average of v(T) over a 
range of T < T c , without invalidating the analytical result at precisely T = T C . 

2.2 Surface fluctuations 

There is another way to consider the spatial fluctuations of hadronization sites 
inside a bin. Imagining the 2D surface to be a membrane that can have dis- 
placements normal to a flat reference surface, then an outward protrusion from 
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the plasma interior can be identified as a bump where hadronization is more 
likely to occur than a dent, an inward indentation. The two competing mech- 
anisms that control the nature of the surface fluctuations are (a) the thermal 
fluctuations, and (b) the confining potential on the partons that prefers no 
deformation of the flat membrane. For the latter we can consider an addi- 
tional term to the free energy, F Sl representing an increase due to spatial 
displacement. Specifically, we take it to be F s = CX)<ij> \zi — Zj\, where C 
is a parameter proportional to the confinement strength, and zt is the vertical 
displacement from a flat surface at site i. C contains the thermal effect rep- 
resented by (/cgT) -1 . The competition between the two mechanisms is then 
introduced by the Boltzmann factor exp(— F s ). 

In our calculation cl we use the GL formalism to describe the mean result, 
which is to be convoluted with the fluctuation result that is computed on a 
lattice by Monte Carlo simulation using the exp(— F s ) factor to determine the 
probability of nonzero Zi displacements. Consideration of cells in bins is again 
necessary, as in the Ising problem described above. Omitting the details we 
summarize here that -F-scaling as in (1) is still valid; moreover, the slopes (3 q 
are nearly independent of C in the range < C < 2, which includes the values 
that has the maximum effect due to surface fluctuations. The formula in (2) 
is also still valid, and the value of v is found ta.be 1.306 ± 0.035. 

Our conclusion is then that the GL resultEm for the scaling exponent v for 
quark-hadron phase transition is unaffected, when spatial fluctuations due to 
surface displacements are taken into account. 



3 Chaotic Behavior 

We now come to fluctuations of a different nature. For phase transitions, 
because of the competition between the ordered (collective) and disordered 
(thermal) motions of the constituents at the critical point, there can be large 
fluctuations over extended domains. But for jets of hadrons produced in hard 
processes or for multiparticle production at low px in hadronic collisions, we 
do not have dense thermal systems of partons (unless it is a heavy-ion colli- 
sion at high energy) and phase transition is not the sort of physics that one 
is concerned with. There are, nevertheless, fluctuations in the event multiplic- 
ity and in the phase-space distribution of the particles produced. Are those 
fluctuations unpredictable? In what way can they be recognized as chaotic 
behavior? 

In order to introduce a quantitative measure of unpredictability we shall 
proceed in stages. First, we shall discuss the notion of erraticity, which con- 
cerns the fluctuations of the spatial patterns of the final states in momentum 
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space. Then we shall discuss entropy indices, followed by their applications 
to both the classical nonlinear problems and the QCD parton showers. The 
consideration involves an area of phenomenology that has hitherto been unex- 
plored. 

3. 1 Erraticity 

Intermittency, refers to the scaling behavior of the multiplicity fluctuations in 
bins of size <5u In determining the normalized factorial moments by averaging 
over all events at a fixed bin and then averaging over all bins, information 
on the spatial patterns from event to event is lost. Erraticity analysis is an 
attempt to capture that information. 

Let F q denote a specific quantification of the spatial pattern in the phase 
space of a final state. It can be the horizontal moments, the correlation in- 
tegrals, or the result of wavelet analysis. After many events one obtains a 
distribution P{F q ) of F q . Now, let us define the normalized (vertical) mo- 
ments! 

Cp, q = (F>)/(F q )* , (6) 

where the averages are done with P(F q ) as the probability distribution. The 
fluctuations of F q become more simple to categorize and more interesting to 
study, if C Pt q exhibits scaling behavior, i.e., 

C M «JV^W , (7) 

where M is the numbeiiof bins in a fixed space. For convenience, we refer to this 
behavior as erraticityc3 It is a natural extension of the notion of intermittency. 
The value of p can be any positive real number. It should not be negative 
because F q may vanish for some events, so P(F q = 0) ^ for some q. For 
Q < p < 1, it is the F q < 1 region that is probed by ip q (p), while for p > 1 
it is the spiky events with high F q that tp q (p) describes. In practice it is not 
necessary to consider p greater than 3. 

Quantities similar to C Pig haiffi been considered before. In statistical 
physics the random energy model liil for spin glass has been investigated with 
the consideration of the quantity (Z p (f3)), where Z(/3) is the partition function 

g w exp [—(3E(uj)] , (3 being the inverse temperature. It led Brax and Peschanski 
to study in the a model the quantity (Z?) where Z q = ^2 m (p m / Ylm P™) 9 > 
p m being the density of particles in the mth bin. Since q plays the role of 
j3 in Z(f3), the possibility of "non-thermal" phase transition is considered in 
Ref. 12. The same quantity (Z^) has similarly been studied in Ref. 13 also 
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in the framework of the a model. The emphasis of erraticity analysis is not 
on the theoretical possibility of a phase transition, but on determining the 
scaling behavior of C p . q from data (experimental and model simulation) and 
on extracting the erraticity indices ip q {j>). It is reasonable to suggest that 
tpqip) provides a stringent test of the reality of any model on multiparticle 
production. 

The study of erraticity can be applied to many problems in physics, chem- 
istry and beyond. Fluctuation in spatial patterns is a ubiquitous phenomenon. 
The Ising model, for instance, has many types of patterns in the lattice space 
depending on the temperature. One can even determine the critical tempera- 
ture for a finite lattice in any dimension by studying the behavior of ipq(p)c-i 

It is also possible to determine the erraticity spectrum e(a) analogous to 
the multifractal spectrum f(a). The definition for e(a), for a fixed q, isll^l 

e(a) = pa — ip(p) , a = dip (p)/ dp . (8) 

At p = 1, e(a) = a, since "0(1) = 0. That value of the spectrum has a particular 
significance, as we shall discuss next. The function e(a) exhibits that value 
explicitly, unlike 



3.2 Entropy indices 

The unpredictability of the outcome of a dynamical process should be related 
to some quantity like the entropy that measures the ignorance about the sys- 
tem. We do not want to follow the time evolution of the system, which is not 
observable. Our entropy must refer to an ensemble of final states. To that end 
we use F q as in Sec. 3.1, to describe the spatial pattern of an event, and P{F q ) 
to describe the distribution of F q after N events. The average of any function 
f(F q ) can also be written as 

M 

dF q P{F q )f{F q )=M- 1 Y,f{F e q ) , (9) 

e=l 

where F* is the F q for the eth event. Introduce now 

M 

P e = FI/Y;F e q , (10) 

e=l 

in terms of which we can define an entropy in the "event space" 

S=-J2 P e ] *Pe , (11) 
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(which should probably be called "eventropy" , since it is not the usual entropy.) 
To calculate S it is more convenient to introduce the moments H p such that 



S = - —lnRp 
dp 



P =i 



From (|9|) we have 
H P =M J dF q P(F q ) 



P = M(Ff) 



(12) 



(13) 



[AfJdF q P(F q )F q _ 
If the system exhibits erraticity, it then follows from (0) that 

5 = m(A/M-^) , (14) 

where 

d 



(15) 



P =i 



We refer to fi q as the entropy indicesta 

If fj, q = 0, then S = hxA/", which is large. One may think that it describes 
a system that is chaotic. Actually, the opposite is true. It corresponds to 
P e = 1/M in jll|), i. e., every event has the same value of F qi according to 
(|l0|). That is hardly what one would expect of a chaotic system, since the 
spatial patterns are the same for every event. S is large because in the event 
space F q is evenly distributed over the entire space. S would be small if F q 
is restricted to only a few events; that would mean large fluctuations in F q . 
Thus finite, nonvanishing positive values of fi q corresponds to wide P(F q ), 
which in turn means unpredictable spatial pattern from event to event. In 
short, fi q is a measure of unpredictability. 

An alternative way of calculating [i q is to circumvent S and work directly 
with C p , 9 by expressing it as 



(F q ) 



Then we have 



dp p ' q 



(16) 



(17) 
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On the other hand, if C PtQ has the scaling behavior (Q), then we also have 

= (j, q In M (18) 

P =i 

in the scaling range of M. Consequently, we obtain 

which is to be determined in the region where a q = ($ 9 In <&„) exhibits a linear 
dependence on InM. 

3.3 Chaoticity 

We now consider the question whether multiparticle production processes are 
chaotic. Let us fixst recall the properties of chaotic behavior in classical non- 
linear dynamicsu-3 Since in such systems trajectories in space-time exist, the 
distance function d(t) between two trajectories can be defined. A system is 
chaotic if d(t) ~ e At , A > 0, no matter how small d(0) = e may be, for e > 0. 
For classical Yang-Mills dynamics the Lyapunov exponent A has been shown 
to be positive by lattice calculationllB 

For quantized Yang-Mills fields the problem is vastly more complicated. In 
addition to the ambiguity associated with quantum chaos in the realm of first 
quantization, we now have also the problems of nonconservation of the number 
of degrees of freedom and of the lack of a meaningful definition of trajectory. 
The absence of an unambiguous notion of time in production processes further 
results in the nonexistence of the Lyapunov exponent for the problem. 

Our first task must be the search for a measure of chaos appropriate to 
our problem. Since nonperturbative QCD is too difficult to implement, we 
focus on parton showers in pQCD, i. e., quark and gluon jets in a tree-diagram 
approximation of QCD branching processes. In place of classical trajectories 
whose initial points are all in a small neighborhood of one another, wc consider 
branching processes all starting from exactly the same virtuality Q 2 and let 
quantum fluctuations take them to different final states. Since we have no 
experimental access to the process of branching, and theoretically the degra- 
dation of vituality does not have an unique association with the temporal 
evolution, our measure of chaos must be focused on what we can observe, viz, 
the characteristics of the final states. It is in that respect that the subject of 
erraticity and entropy indices becomes relevant. 

There are two issues here. One is to find a measure of chaos, and the other 
concerns QCD branching processes. We discuss them separately. 



dp p ' q 
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3.3.1 A new measure of chaos 

Our proposal is to use fi q as an alternative measure of chaos in problems 
where only the spatial patterns can be observed. Whether that approach 
agrees with the conventional approach can be determined only by applying 
the measure to knopii chaotic systems, of which there are many with simple 
nonlinear dynamicsEj They are classical systems with well-defined trajectories, 
and the positivity of the maximum Lyapunov exponent A is a sufficient criterion 
for chaos. Since our measure is concerned with spatial patterns rather than 
temporal evolution, it is necessary to construct a spatial pattern for every 
trajectory. That is achieved by considering a set of points corresponding to 
the positions of a trajectory at a discrete collection of times, separated by finite 
intervals apart. Clearly, two nearby trajectories for A < would generate two 
very similar patterns, while two chaotic trajectories for A > would have very 
dissimilar patterns. We can generate M events by considering M trajectories 
all starting from a small neighborhood of an initial point. Studying the P{F q ) 
distribution of those events and then comparing the resultant fi q to A constitute 
our procedure to check the effectiveness of [i q as a measure of chaos. 

We have applied this procedure to the logistic map and the Lorenz at- 
tractor, the details of which are omitted here. They are nonlinear systems 
that become chaotic when their control parameters r exceecLcertain critical 
valuesll3 We present here only the result for the logistic mapEJ In Fig. 1 the 
dashed line shows the value of A as a function of r. For r > 3.57, A becomes 
positive, although there are short intervals where A drops below zero. The solid 
line indicates the value of (12 (normalized by a specific factor), calculated at 
discrete points in the same range of r. Evidently 112 coincides very well with A 
except that it cannot be negative but vanishes when A < 0. The multiplicative 
factor used for /12 in the plot depends on the number of spatial points taken 
to determine the P(F q ) distribution, and has no particular significance. What 
is of great significance is that we now have an alternative measure of chaos. 
The positivity of A is coordinated with the positivity of ^2 (and all other fi q ), 
so that the fluctuations of spatial patterns can equally serve as a means for 
revealing chaotic behaviors. 

There are many complex systems possessing complex patterns. The de- 
termination of their entropy indices may reveal certain features that are not 
otherwise recognizable. For some nonlinear systems that are known theoret- 
ically to be chaotic, the experimental verification of A may be difficult, since 
the precise adjustment of the initial condition may not always be possible. 
Our method of studying the spatial patterns may therefore offer a feasible 
alternative. 
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3.3.2 Chaos in QCD parton showers 



Since perturbative QCD is well developed theoretically, and the final states of 
quark and gluon jets are well measured experimentally, we propose that the 
erraticity of hard processes be used as a new arena for comprehensive tests 
and that the entropy indices of parton showers be determined to reveal the 
chaoticity of QCD. 

Since hadronization is irrelevant to the question of chaos in pQCD, we have 
developed an efficient Monte Carlo generator of quark and gluon jets with the 
usual Sudakov form factor and splitting functions incorporated in the codec! 
The initial virtuality of a quark (Q jet) or a gluon (G jet) is fixed at a Q 2 
value for all events, and all partons evolve and branch until the virtuality of 
a parton reaches < Qg, where branching is terminated for that_parton. From 



which the inclusive distribution is flat), we calculate the horizontal normalized 
factorial moments F q , then the distribution P(F q ), the vertical moments C Ptq , 
the erraticity indices ip q (p), and finally the entropy indices fi q . The results of 
the calculation for Q/Qo = 10 3 are shown in Fig. 2. The positivity of pi q for 
both Q and G jets indicates that the QCD branching processes are chaotic. 
Moreover, the Q jets are more chaotic than the G jets. The reason is that quark 
jets have fewer partons produced, so the multiplicity fluctuations relative to 
the mean is larger. 

We have also considered the fixed coupling problem, since by varying a s 
as a control parameter one can study the behavior of the system at the onset 
of chaos. It turns out that there is no threshold of chaos. When a s is small, 



the spatial pattern of each event 
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the multiplicity is so small that fi q is large for both Q and G jets. Thus 
chaoticity is an inherent property of QCD dynamics and cannot be tuned out 
by decreasing a s , the only tunable parameter that we can associate directly 
with the strength of the nonlinear term. 

Using the positivity of [i q as a criterion for chaos, we have not only found 
that the non-Abelian, nonlinear dynamics of QCD is chaotic, but also come to 
realize that any quantum system in which the number of quanta of interest is 
not conserved is likely to have positive \x q . Thus while chaotic behaviors are 
spectacular and remarkable phenomena in classical dynamics, their occurrences 
seem to be common and generic in quantum systems. Perhaps the notion of 
chaos should not be more emphasized in multiparticle production than the 
realization that the unpredictability of the final states can be quantified. To 
calculate the entropy indices for all collision processes, soft as well as hard, 
that can agree with the experimental data would be a theoretical challenge. 
Our results obtained should not be taken quantitatively for comparison with 
experiments, since the MC code has not been tuned to check the other features 
of parton showers. However, they have been effective in demonstrating the 
procedure of determining fi q and in elucidating the question of chaos in QCD. 

3.4 Hadronic and nuclear collisions 

For hadronic collisions we expect more fluctuations than in e + e~ annihilation, 
since there is impact-parameter variation from event to event. Even at fixed 
impact parameter there can be various cut-Pomerons, each of which can have 
large variations in the final states corresponding to many possible ways that 
the branching processes can proceed, at least according to the geometrical 
branching model (GBM) for soft processeslia That does not even include the 
consideration of hard subprocesses at high energies that would further increase 
the degree of fluctuations. Gianini has found in_his preliminary analysis of the 
old Fepnilab data that is in excess of 0.4t£l Cao is currently upgrading 
ECCOpI an event generator that is based on GBM, by taking into account 
resonance production. The preliminary result on the entropy indices is that 
they are also very large. 

In heavy-ion collisions the fluctuations in impact parameter are usually 
controlled by making cuts in Ex- Because of the high multiplicity per event 
not much has been observed in the intermittency behavior. More may be 
revealed in the erraticity analysis. I venture to speculate here what the results 
on \x q would be in the four possible scenarios: with and without a hadron-gas 
phase, and with and without a quark-to-hadron phase transition. If there is a 
hadron gas phase before hadrons are emitted so that there is thermalization 
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in the final state, then the emitted hadrons are randomized. Random spatial 
patterns of many particles can differ significantly from event to event, so I 
would not expect fj, q to be small whether or not there is a phase transition. If 
there is no hadron gas phase that randomizes the hadronic final states, then one 
should expect to see the dynamical effects of the production processes in the 
erraticity analysis. If there is a phase transition, the produced hadrons should 
exhibit long-range correlations that result in clusters of all sizes (assuming 
second-order phase transition and Kadanoff scaling). In that case we should 
expect larger values of \i q than in the case of no phase transition. 

4 Conclusion 

Significant progress has been made recently on the subject of fluctuations in 
multiparticle production. On critical behavior considerations of fluctuations 
beyond the result of mean-field theory have not altered the scaling exponent v. 
For noncritical phenomena new measures of fluctuations from event to event 
are proposed. Experimental determination of v, ^P q {p) and fj, q is urgently 
needed. For v it is necessary to consider high q moments. For tp q (j>) and [i q 
low q values like 2 and 3 are enough and p in the range 0.5 to 2 is sufficient. 
It is highly likely that many existing event simulators may not be able to 
reproduce the experimental data on those quantities. Significant change of p q 
due to quark-gluon plasma formation would be very exciting. 
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